Local and Global Casimir Energies for a Semitransparent 

Cylindrical Shell 

Ines Cavero-Pelae^] and Kimball A. MiltorQ 

Oklahoma Center for High-Energy Physics and Homer 
L. Dodge Department of Physics and Astronomy, 
University of Oklahoma, Norman, OK 73019-2061 

Klaus Kirsteif] 

O Department of Mathematics, Baylor University, Waco, TX 76798 

^ ! (Dated: February 1, 2008) 

Oh- Abstract 
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! ated with a A<5-function potential in a 3+1 dimensional circular cylindrical geometry are considered. 
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cal energy density only when the latter is supplemented by an energy term residing precisely on 
the surface of the cylinder. The latter is identified as the integrated local energy density of the 
cylindrical shell when the latter is physically expanded to have finite thickness. Inside and outside 
the 5-function shell, the local energy density diverges as the surface of the shell is approached; the 
divergence is weakest when the conformal stress tensor is used to define the energy density. A real 
global divergence first occurs in 0(A 3 ), as anticipated, but the proof is supplied here for the first 
time; this divergence is entirely associated with the surface energy, and does not reflect divergences 
■ in the local energy density as the surface is approached. 
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I. INTRODUCTION 



The subject of self-energies due to quantum fluctuations of fields constrained by physical 
boundaries, such as the Casimir energy due to a perfectly conducting spherical shell has 
been controversial nearly from the outset 0, HI . Partly this is because such a self-energy is 
apparently not well-defined, for it is not accessible by deforming part of the surface. More 
fundamentally, it is because there are strong divergences in the local energy density as the 
(idealized) boundary is approached, which would seem to rule out the existence of a finite 
total energy for the system. How are such divergences to be squared with the mathematical 
proof that the electromagnetic Casimir energy of a region with a closed, smooth, perfectly 
conducting boundary is finite 0? 

These issues have been brought to the forefront by a series of papers by Graham et 
al. [Hf. Essentially, they assert that it is impossible to ascribe any physical meaning to the 
self-Casimir energy of a single object, such as Boyer's sphere 0. However, the divergence 
issues they raise are hardly new for example, for a massless scalar particle in the presence 
of a spherical 5-shell potential, governed by the Lagrangian 



C = --d^d^ - --8(r - a)0 2 , (1.1) 



2 ^ y 2 a 2 

a divergence occurs in third-order in A 0, and possible ways of dealing with it have been 
suggested [1, E|. 

For a different viewpoint see the considerations of Barton [10(. He argues, that for 
any single connected body the physics is dominated by the divergent components of the 
self-energy. The pure Casimir terms, which in the examples studied in Ref. jjjj] are the 
convergent components when the no-cutoff limit is taken, are shown to be much smaller 
such that they will never become observable. 

In a recent paper E3 the model described by (jl.lj) was considered in some detail. In 
fact, the model examined was somewhat more general, in that the 5 function was replaced 
by a step-function potential of width S and height h. As 5 — > and h — > oo with hS held 
fixed at unity, a 5-function potential is recovered. For such a potential we calculated the 
energy density inside and outside the region of the potential, from which we could calculate 
the total energy as a function of A, as well as the local energy density. As long as the 
potential is finite, the local energy density may be integrated to yield the total energy. In 
the singular 5-function limit, however, the energy within the shell becomes a localized surface 
energy which must be added to the integrated local energy of the regions inside (r < a) and 
outside (r > a) the sphere, which has precisely the anticipated form [l^, EH EH Ell- For 
weak coupling with the 5-shell potential, the total Casimir energy is finite in 0(X 2 ), but 
divergent in third order, which divergence precisely corresponds to the divergence of the 
surface energy in that order. Thus, it is plausible that such a divergence should be absorbed 
in a renormalization of the surface energy. 

In this paper, we turn to the corresponding cylindrical case. The situation in many ways 
is similar. However, there are curiosities associated with the cylindrical geometry that make 
this new analysis intriguing. First, the Casimir self-stress on a perfectly conducting circular 
cylinder was found [16[ to be attractive, and of somewhat smaller magnitude, compared 
to the repulsive stress found by Boyer for a perfectly conducting sphere hi. It was found 
that a dilute dielectric cylinder had vanishing van der Waals energy [It], Il8l |. as did a dilute 
dielectric-diamagnetic cylinder (with the speed of light the same on the inside and the 
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Type -EsphcrcQ £c y iindcra 2 References 



EM 


+0.04618 


-0.01356 


m [is] 


D 


+0.002817 


+0.0006148 


[30j[31] 


-I) 2 


+0.004767 = , 23 fi 





[25j[21] 


e 


+0.04974 = 4- 

OZ7T 





[22] [ia 


be 2 


±0.0009 





[29] [23] 


A 2 


+0.009947 = 

OZ7T 





[32] 



TABLE I: Casimir energy (E) for a sphere and Casimir energy per unit length {£) for a cylinder, 
both of radius a. Here the different boundary conditions are perfectly conducting for electro- 
magnetic fields (EM), Dirichlet for scalar fields (D), dilute dielectric for electromagnetic fields 
[coefficient of (e — l) 2 ], dilute dielectric for electromagnetic fields with media having the same 
speed of light (coefficient of £ 2 = [(e — l)/(e + l)] 2 ), for perfectly conducting surface with eccen- 
tricity 5e (coefficient of <5e 2 ), and weak coupling for scalar field with 5- function boundary given by 
Ql.l|) . (coefficient of A 2 ). The references given are, to the authors' knowledge, the first paper in 
which the results in the various cases were found. 



outside) HI H 0, which seemed to imply the vanishing of the Casimir ener gy f or a 
dielectric cylinder of permittivity e in order (e — l) 2 , as was only recently verified |2l|, l22j . 
If a perfectly conducting cylinder is slightly deformed by giving its cross section a slight 
eccentricity 8e, the change in the Casimir energy vanishes in order 5e 2 (HJ. None of these 



vanishings occur for a sphere |24J, |2jj 120, |27J, UB, l2i|- so they all reflect the flatness of 



the cylindrical geometry. In this paper we establish another example of the second-order 
vanishing effect for the cylinder, that is, that the 0(X 2 ) term in the Casimir energy for the 
semitransparent cylinder is zero. A summary of the facts comparing results for sphere and 
cylinder, together with the first reference for each result, is given in Table |U 

The outline of this paper is as follows. In the next section we will derive the Green's 
function for the semitransparent cylinder. In Sec. IIHI we will then compute the Casimir 
pressure on the cylinder, and thereby infer the Casimir energy. In Sec. IIVI we will rederive 
that energy directly. The weak-coupling evaluation of the Casimir energy will be the subject 
of Sec. |3 We will show that the Casimir energy through order A 2 vanishes both by analytic 
and numerical arguments, with an explicit isolation of the divergent term which may be 
unambiguously removed. An independent derivation of these results using zeta function 
techniques is given in Sec. IVII The proof that for a cylinder of arbitrary cross section the 
divergence in the Casimir energy occurs in 0(X 3 ) is supplied in Sec. IVHI As we have noted in 
Table HI the strong-coupling (Dirichlet) result was earlier derived by Gosdzinsky and Romeo 
[3ll | ; we reproduce their result in Sec. IVII II Again the explicit divergent terms that must be 
removed by renormalization are identified. We turn to an examination of the local energy 
density in Sec. 11X1 We show that the integrated local energy density differs from the total 
energy by a surface term that resides precisely on the cylindrical surface. We also examine 
the surface divergences that appear in the energy density as one approaches the surface; the 
leading divergence, which is independent of the shape of the surface, may be eliminated by 
choosing the conformal stress tensor. Then one is left with surface divergences, in either 
strong or weak coupling, that are exactly one-half that for a sphere, reflecting the vanishing 
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of one of the principal curvatures for a cylinder. In Sec. [Xj we examine what happens when 
the 5-function potential is thickened to a cylindrical annulus of thickness S and height h. 
We show that then the surface energy is resolved as the integrated local energy density of 
the field within the confines of the annulus. The divergence in the total energy in 0(X 3 ) 
for a 5-function shell is exactly that due to the surface energy alone. Finally, we offer 
a perspective of the situation in the Conclusion. Appendix A offers further details on the 
divergences occuring in the global theory, particularly in 0(A), while Appendix B elaborates 
some further aspects of the surface energy. 

II. GREEN'S FUNCTION 

We consider a massless scalar field in a (^-cylinder background, 

Ant = ~5(r - a)0 2 , (2.1) 

a being the radius of the "semitransparent" cylinder. We recall that the massive case was 
earlier considered by Scandurra |33| • Note that with this definition, A is dimensionless. The 
time-Fourier transform of the Green's function, 

G(x,x') = j |V^'>£(r,r'), (2.2) 

satisfies 

£(r,r') = <Kr-r'). (2.3) 
Adopting cylindrical coordinates, we write 

£(r,r')= / ^e^" 2 ') V ^-e^^g m (r, r'; k), (2.4) 

m=— oo 

where the reduced Green's function satisfies 

9m(r, r'; k) = -5(r — r'), (2.5) 
r 

where k 2 = k 2 — uo 2 . Let us immediately make a Euclidean rotation, 

uj -> <, (2.6) 

where ( is real, so k is likewise always real. Apart from the 6 functions, this is the modified 
Bessel equation. 

Because of the Wronskian satisfied by the modified Bessel functions, 

K m {x)I' m {x) - K' m (x)I m (x) = (2.7) 

x 

we have the general solution to ()2.5|) as long as r ^ a to be 

g m (r, r'; k) = J m (/tr < )K m (/€r > ) + A(r')I m ( K r) + B(r')K m ( K r), (2.8) 



-V 2 -co 2 + -S(r - a) 
a 



1 d d 2 m A . . 

r r ~r + K H — t + _ "( r _ a ) 

r dr dr r z a 
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where A and B are arbitrary functions of r' . Now we incorporate the effect of the 8 function 
at r = a in (|2.5|) . It implies that g m must be continuous at r = a, while it has a discontinuous 
derivative, 

r=a+ 

= \g m {a,r';k), (2.9) 



d 



a fo9m(r, r'; fc) 

from which we rather immediately deduce the form of the Green's function inside and outside 
the cylinder: 



r,r'<a: g m (r,r';k) = I^kv^K^kv^ 



— — — — - -I m {nr)L m {nr ), 

1 + Al m {K,a)K m {K,a) 



(2.10a) 



r,r'>a: g m (r,r';k) = I m (/tr < )i ; C m (Kr > ) 

All Oca) 



1 + XI m (Ka)K m (m) 



K m (Kr)K m (Kr'). (2.10b) 



Notice that in the limit A — > oo we recover the Dirichlet cylinder result, that is, that g n 
vanishes at r = a. 



III. PRESSURE 

The easiest way to calculate the total energy is to compute the pressure on the cylindrical 
walls due to the quantum fluctuations in the field. This may be computed, at the one-loop 
level, from the vacuum expectation value of the stress tensor, 



d^d' u - ^ v d x &^j ~G(x,x') 



£{d»d u -g» v d 2 )-G{x,x). (3.1) 



Here we have included the conformal parameter £, which is equal to 1/6 for the stress tensor 
that makes conformal invariance manifest. The conformal term does not contribute to the 
radial-radial component of the stress tensor, however, because then only transverse and time 
derivatives act on G(x,x), which depends only on r. The discontinuity of the expectation 
value of the radial-radial component of the stress tensor is the pressure of the cylindrical 
wall: 

P (-^rr/in (-^rr/out 

Xk 2 



-| v /"OO poo 



167r3 m _ 0Q J-oo J-oo 1 + \I m (Ka)K m (Ka) 
x [Kl(Ka)I' 2 (Ka) - ll(Ka)K' 2 (Ka) 



1 00 poo poo J 

— J dk J d(-—ln[l + \I rn (Ka)K rn (m)}, (3.2) 



m=— oo 



167T 3 

where we have again used the Wronskian (|2.7J) . Regarding ka and (a as the two Cartesian 
components of a two-dimensional vector, with magnitude x = Ka = a 2 , we get 

the stress on the cylinder per unit length to be 

S = 2naP = / dxx 2 > — In [1 + \I m (x)K m (x)] , (3.3) 

Ana 6 J Q / — 4 dx 



m=—oo 
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which possesses the expected Dirichlet limit as A — > oo. The corresponding expression for 
the total Casimir energy per unit length follows by integrating 



that is, 



I roc 00 d 

£ = -!-^ dxx 2 V — \n[l + XI m (x)K m {x)]. (3.5) 
8na z Jr. ax 

m=—oo 



This expression is of course, completely formal, and will be regulated in various ways in the 
following, for example, with an exponential regulator (in Sees. W \ IVIII[ Appendix A), or by 
using zeta-function regularization (in Sec. IVI|) . 

IV. ENERGY 

Alternatively, we may compute the energy directly from the general formula 0] 

To evaluate (|4.1j) in this case, we need the indefinite integrals 

J*dyyl 2 m (y) = \[{x 2 + m 2 )ll{x)-x 2 C], (4.2a) 

r°° i 

J dyyK 2 m (y) = --[(x 2 + m 2 )K 2 m (x)-x 2 K' 2 }. (4.2b) 

When we insert the above construction of the Green's function, and perform the integrals 
as indicated over the regions interior and exterior to the cylinder we obtain 

2 °° poo poo i j 

£ = -7T2 E / d W dk( 2 -—\n[l + \I m (x)K m (x)]- (4.3) 

o7T _ J —oo J —oo 



m=—oo 



Again we regard the two integrals as over Cartesian coordinates, and replace the integral 
measure by 



OO POO POO 

2 _ „ / J„ „3 



d( dk( 2 = 7i dn k . (4.4) 

oo J —oo J 



The result (13.5)1 immediately follows. 



V. WEAK-COUPLING EVALUATION 



Suppose we regard A as a small parameter, so let us expand ()3.5|) in powers of A. The 
first term is 

\ 00 f°° d 
£ {1) = -tt^ E / dxx 2 ^K m {x)I m {x). (5.1) 
8na 2 Jr. dx 

m=—oo w 
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The addition theorem for the modified Bessel functions is 

oo 

K (kP)= e m ^K m (kp)I m (kp), P > p' , (5.2) 

m=— oo 

where P = a/p 2 + p' 2 — 2pp' cos(4> — <f>') ■ If this is extrapolated to the limit p' = p we 
conclude that the sum of the Bessel functions appearing in (j5.1j) is K (0), that is, a constant, 
so there is no first-order contribution to the energy. For a rigorous derivation of this result, 
see Sec. IVI Al and also Appendix A. 



A. Analytic Regularization 

We can proceed the same way to evaluate the second-order contribution, 

£ {2) = 77T- 2 / d^ 2 ^ E ( 5 - 3 ) 
lo7rGr In ax z — ' 

u m=— oo 

By squaring the sum rule (|5.2jl . and taking the limit // — > p, we evaluate the sum over Bessel 
functions appearing here as 

I 2 m (x)K 2 m (x)= £K 2 (2xsmcp/2). (5.4) 



m=— 00 



Then changing the order of integration, the second-order energy can be written as 



where the Bessel-function integral has the value 1/2. However, the integral over (p is diver- 
gent. We interpret this integral by adopting an analytic regularization based on the integral 



21 



2t „ 

if \ s 



which holds for Res > —1. Taking the right-side of this equation to define the <p integral 
for all s, we conclude that the ip integral in (|5.5j) . and hence the second-order energy 8^ 2 \ 
is zero. 



B. Numerical Evaluation 

Given that the above argument evidently formally omits divergent terms, it may be more 
satisfactory, as in [2]]], to offer a numerical evaluation of £^ 2 \ (The corresponding argument 
for is given in Appendix A.) We can very efficiently do so using the uniform asymptotic 
expansions (m — > 00): 



^-(i+D-i)'^). (™) 
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where x = mz, t = 1 / y/1 + z 2 , and the value of r\ is irrelevant here. The polynomials in t 
appearing in (|5.7|) are generated by 



u {t) = 1, u k {t) = -t 2 (l-t 2 )u' k _ 1 (t) + - ds(l-5s 2 )u k ^(s). 



(5.8) 



Thus the asymptotic behavior of the product of Bessel functions appearing in (J5.3J) is 



I m {x)K m {x) 



Am? 



r k {t) 



k=l 



m 



2k 



The first three polynomials occurring here are 



n{t) 



t z 



16 
f_ 
16 



;i-et 2 + 5t 4 ), 

(7 - 148t 2 + 554t 4 - 708t 6 + 295t 8 ), 
(36 - 1666t 2 + 13775t 4 - 44272t 6 



+ 67162t 8 - 48510t 10 + 13475t 
We now write the second-order energy ()5.3j) as 



12\ 



x- 



dx x 



ll{x)Kl(x] 



1 



l v /l f \ 
- hm - + > m 
4 s^o \ 2 ^ ) 

POD 



4(1 + x 2 ) 
d 1 



dz Z 2 ~ S - 

I, dz 1 + 



2/V 



m=l fc=l 



m=l 



£(s)j£(aO - 



£ 2 



4m 2 



r k (t) 



k=l 



III 



2k 



In the final integral z = x/m. The successive terms are evaluated as 

A 2 



£(2) 



1/ , ,n 1, n C(2 7( 4 31C 6 

-(7 + In 4) - -ln2?r - + 

4 W ; 4 48 1920 16128 



87ra 2 

+0.000864 + 0.000006 



A 2 



■(0.000000), 



(5.9) 



(5.10a) 
(5.10b) 

(5.10c) 



(5.11) 



(5.12) 



where in the last term in (|5.11j) only the m = 1 and 2 terms are significant. Therefore, we 
have demonstrated numerically that the energy in order A 2 is zero to an accuracy of better 
than 10~ 6 . 

The astute reader will note that we used a standard, but possibly questionable, analytic 
regularization in defining the second term in (|5.11|) . where the initial sum and integral are 
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only denned for 1 < s < 2, and then the result is continued to s = 0. Alternatively, we could 
follow Ref. 2l( and insert there an exponential regulator in each integral of e~ xS , with 5 to 



be taken to zero at the end of the calculation. For m/fli: becomes mz, and then the sum 
on m becomes 

oo 

Y, e ~ mzS = ^Y- ( 5 - 13 ) 



m=l 



Then when we carry out the integral over z we obtain for that term 

£-ilnar. (5.M) 
Thus we obtain the same finite part as above, but in addition an explicitly divergent term 

£ - = (5 - 15) 

If we think of the cutoff in terms of a vanishing proper time r, 5 = r/a, this divergent term 
is proportional to 1/a, so the divergence in the energy goes like L/a, if L is the (very large) 
length of the cylinder. This is of the form of the shape divergence encountered in Ref. jijj . 



VI. ZETA-FUNCTION APPROACH 

For the massless case, in the zeta-function scheme the regularized total energy reads 

£ = --— a 2s (l + *[-! + 2 ln(2/x)]) V / dxx 2 - 2s — ln(l + XK m (x)I m (x)) , (6.1) 

oTTCL Jq OLX 



where /i is an arbitrary mass scale. This result can, for example, be taken from Ref. [33]; 
see also (8.3.16) of Ref. j35|. In order to find the total energy one needs to find the analytic 
continuation of this expression to s = 0. Formally, in the s — > limit, this is exactly (J3.5j) . 
so this is the zeta-function regular ization of that result. 

A. 

The first order of the energy, in the weak-coupling expansion, reads 

A °° f°° H 

£ {1) = -—^a 2s (l + s[-l + 2\n(2^]) V / dx x 2 ' 2s —K m {x)I m {x) 
6ira z ^-^ In ax 

m=—oo u 

= --^a 2s (l + S [-l + 21n(2/x)]) / dx x 2 ~ 2s —K (x)I (x) 
87ia z I J dx 

°° f°° d 1 

+ 2^2m 2 ~ 2s dz z 2 - 2s —K m {mz)I m {mz)y (6.2) 



rn=l 



where, for the purpose of the following argument, the m = contribution has been separated 
off. To be precise, let us mention that the m = contribution should be thought of as coming 
from a limit M — > of the integral . This is what one actually obtains in the contour 
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integral formalism; for full details see Ref. 35||. With this in mind, (6.2) is analytic in the 
strip 1 < Res < 3/2. The substitution x = mz has been done in the last integral to make 
the convergence properties of the series and the integral better apparent when using the 
uniform asymptotics of the Bessel functions (as noted above; see below). 

We want to find the analytic continuation of this expression to s = 0. As it stands, the 
expression is not valid about s = and manipulations need to be done. In order to evaluate 
it further, one would like to interchange summation and integration and use the addition 
theorem for the Bessel functions. From the uniform asymptotic expansion of the Bessel 
functions (m > 0), 



K m (mz)I m (mz) 



~— + ^(l-6t 2 + 5t 4 ) + G 
2m 16m d 




(6.3) 



as we have seen already in ()5.9|) and (j5.10a|) . it is clear that the resulting summation is 
divergent at s = 0; one can see that as it stands the summation is convergent only for 
Res > 1. One therefore needs to subtract and add as many asymptotic terms as needed 
to make the resulting summation convergent. In this procedure the hope is always that the 
asymptotic terms which potentially could contain singular contributions at s = can be 
handled analytically. Subtracting all asymptotic terms given above, one sees the resulting 
summation is well defined at s = 0. 

Let us define £^ 1 ' without the prefactor in ()6.2|) . 



go = £ / 



dx x 2 2s -^-K m (x)I m (x), 



dx 



(6.4) 



m=— oo 



which, given the above remarks, we rewrite as 



£{1) 



x d 

2-2s u 



dx x 

oo 

m=l 

oo 

+ 2£ 



dx 



Kn(x)Ux) 



m 



-OO J 

2-2s / j 2-2s a 



711 



2-2s 



rn=l 



dz z l ~ ls ^- 

dz 



00 d 

dz z 2 ~ 2s 4- 
dz 



5t 4 



t t 3 (l — 6t 2 -f 
2m 16m 3 



■ t t\\ - Qt 2 + 5t 4 

K m (mz)I m {mz) - 

2m Ibm 4 



(6.5) 



We need to find the analytical continuation of this expression to s = 0. For the first term 

r(2 - s)r(i - s)v (s - \) 



we find for ~ < Re s < 1 



d 

2-2s u 



dec oo 



dx 



K (x)I (x) 



2v^r(s) 



(6.6) 



where the right side vanishes when analytically continued to s = 0. 

For the middle, asymptotic, terms in ()6.5j) we use for 1 — |<Res<2 



x d 

2-2s u 



dz z 



1 



dz (1 + Z 2 )™/2 



r(2-s)r(f + s-i) 



(6.7) 
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for n an integer. The m-series then leads to the zeta function of Riemann, and putting these 
remarks together we first obtain 



dz 



t t 3 (l - Qt 2 + 5t 4 
2m 16m 3 



<r(2s 



r(2- 8 )r( 8 -|) 



7T 



i , s , f r(s + i) r(s + |) r( s + §) 

+ gfad + 2 S )r ( 2 - .) |_L^ + - 

= - Ch(2s - D r(2 -^> - &(! + 2s) f(iziH|zi)£(f±i) . (6 . 8) 

V 71 " ^V 71 " 
This derivation, which assumes 1 < Res < 2, provides the analytic continuation to s — 0. 
In the second term one needs to be a little careful as the pole at s = coming from 
the Cr(1 + 2s) multiplies s, contributing altogether something finite. Adding up the two 
contributions, namely —1/6 and +1/6, the total contribution is 0. 

This shows that m = and the asymptotic terms do not contribute to order A. The 
fact that the formulas (|6.6j) and ()6.8)1 have non-overlapping domains of validity is irrelevant, 
taking into account the remark made below (|6.2jl . 

Let us now consider the most difficult term (which usually one is not able to handle 
analytically, but here it works out fine). We want to find 

„ m f°° , 2 d A f r , , w , v t t 3 (l - 6t 2 + 5t 4 ) I , . 

which is the last term in ()6.5|) at s = with the substitution x = zm and therefore here 



t = l/yl + (x/m) 2 ; this representation, at least for the first term, is better for the fol- 
lowing application of the addition theorem of the Bessel functions. Remember, we are able 
to put s = and interchange summation and integration because is well defined by 
construction. 

We will now show that this finite quantity actually vanishes. 

As it stands, the summation cannot be performed. In order to be able to deal with the 
individual terms separately, we introduce an oscillatory factor. For Z^ l > we therefore write 

Z« = lim l°° dx x 2 ±2Re V (^ m (x)/ m (x) - ±- - ^ ~ f 3 + ^ ) ■ (6.10) 
v-*°Jn dx ^-^ 2m 16m 3 

m=l v ' 

The advantage is that we can now use the addition theorem for Bessel functions (15.2)1 to see 
that 

oo 

2Re ^e m ^ m (x)J m (x) = K (2xsin |) - K (z)I (z). (6.11) 

m=l 

This shows that 

= lim / dx x 2 — { K ( 2xsin - ) - K (x)I (x) 
¥>->o J dx [ V 2 / 

y oos(m y ) t _ 1 ^ cos(rM _ ^ j 
^ m 8 ^ m 3 v ; J v ; 

m=l m=l J 
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As they stand, the two series cannot be performed further, because t contains a nontrivial 
m-dependence. Of course one would like to substitute back x = zm, but in order to do so 
one needs to separate all terms and perform the x-integration for each term individually. 
This gives divergent results. So in order to be able to consider each term individually, 
we have additionally to regularize the x-integration and will therefore consider (effectively 
reinserting the zeta-function regularization) 



ZQUs) = lim lim [Z u (a, if) + Z 12 (a, if) + Z 13 (a, if) + Z u (a, cp)) , (6.13) 

ip— >0 a^O 



where 



Z n (a,<p) 



[ dx x 2 ~ a -^K (2xsin^y (6.14a) 
Jo dx V 2/ 

f°° d 

Z l2 (a,ip) = - / dx x 2 - a —K (x)I (x), (6.14b) 
Jo dx 

(6.14c) 

J ° m=l 

Z u (a, <p) = -- !°° x 2 ~ a 4- T CQS( T^ 3 (1 - 6t 2 + 5t 4 ). (6.14d) 
8 In dx z — ' mr 

m=l 

The integral for Z\\ can be performed for Re a < 2 and one obtains 

toH -Ir(l-f)r(2-f)sin-(f). (6.15) 



Z„(6, V )) = -ism- 2 ('^ > ). (6.16) 



At a = this gives 

Ziifn.£o) = - . ... , 

2 V2, 

Ultimately we want to send <£> — > 0; remembering that the oscillating factor was introduced 
to make each separate sum convergent, it is expected that in this limit divergences might 
occur. By construction, this divergent piece as (p — > has to be cancelled by one of the 
remaining contributions, since at the beginning the combination of all terms was finite. 
The second contribution is calculated as in (J6.6|) . namely for 1 < Re a < 2, 

Z 12 (a, <p) = (2 - a) [ *> \ 2 j , (6.17) 
which vanishes at a = 0. So this does not contribute to the energy. 



Next we find (we substitute back to the case where t = 1/ a/1 + z 2 ) for 1 < Re a < 4 

Z 13 {a,(p) = -^ dzz 2 - a ^?j ^m^costm^j , (6.18) 



which at a = gives 

which exactly cancels ()6.16|) . 



Zi 3 (0,y)= . s , (6.19) 
2 sin yj/2 
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Finally, we find for — 1 < Re a < 4 



Z M („ „ = -i (f * ^ (1 _ «, + «.)) ^-1) . (6 , 0) 

Looking at the calculation of the asymptotic terms, it is seen that the integral is the same 
as that encountered in the second term in ()6.8|) and hence at a = this vanishes. Note, 
that this time the is not multiplied by any infinite quantity because ip ^ at this stage. 
So no contribution results from this term. 

Adding up all the contributions, the answer is simply £W = and therefore 

£ w = 0, (6.21) 

in agreement with the argument given at the beginning of Sec. |V] (See also Appendix A.) 



B. £( 2 ) 

The second order contribution, fortunately, is significantly simpler. The zeta-regularized 
version of the second-order energy reads 

2s °° r°° A 



,2 



^ poo j 

l + s[-l+21n(2 M )]) £ J dxx 2 -* s 2-Kl(x)P m (x). (6.22) 



167ra 

m=— oo 



As before, let us define the second order term without the prefactor, 

°° POO J 

£ {2) = £ / dxx^f x Kl{x)P m {x). (6.23) 

m=— oo 

As it stands, the representation is well defined in the strip 1/2 < Res < 1. In that strip we 
interchange integral and summation using (|5.4|) . Substituting u = 2xsin<^/2, we obtain for 

EH) 

In the relevant range, namely 1/2 < Res < 1, we have the integral [ ()5.6|) 



o 



#(sin|J = V r ^ ^ - (6.25) 



So the analytic continuation of the integral to s = vanishes. Given that the it-integral is 
well behaved about s = (equaling —1 there), we have already found the final answer to be 
£~W = o and thus 

£ (2) = 0. (6.26) 
The first two orders of the weak-coupling expansion indeed do identically vanish. 
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VII. DIVERGENCES IN THE TOTAL ENERGY 



In the following we are going to use heat-kernel knowledge to determine the divergence 
structure in the total energy. We consider a general cylinder of the type C = K. x Y, where 
Y is an arbitrary smooth two dimensional region rather than merely being the disc. As 
a metric we have ds 2 = dz 2 + dY 2 from which we obtain that the zeta function (density) 
associated with the Laplacian on C is (Res > 3/2) 



Xy Xy 



^(-iW-i'l- (T.i) 



2tt r(s) V 2 

Here Ay are the eigenvalues of the Laplacian on Y, and (y{s) is the zeta function associated 
with these eigenvalues. In the zeta-function scheme the Casimir energy is defined as 



1 



(7-2) 

s=0 



which, in the present setting, turns into 



E Cas = V s 7 %t (y(s - 1) 

2 P 20FT(s-|) S V ' 

Expanding this expression about s = 0, one obtains 



(7.3) 

s=0 



E Cas = ^-Cy(-l) + ^ (Cr(-l) [2 ln(2^) - 1] + &(-l)) + 0(a). (7.4) 

S7TS 87T 

The contribution associated with £y(— 1) can be determined solely from the heat-kernel 
coefficient knowledge, namely 

Cy(-l) = -04, (7.5) 

in terms of the standard 4th heat-kernel coefficient. The contribution coming from Cy-( — 1) 
can in general not be determined. But as we see, at least the divergent term can be deter- 
mined entirely by the heat-kernel coefficient. 

The situation considered in the Casimir energy calculation is a 5-function shell along 
some smooth line E in the plane (here, a circle of radius a). The manifolds considered 
are the cylinder created by the region inside of the line, and the region outside of the line; 
from the results the contribution from free Minkowski space has to be subtracted to avoid 
trivial volume divergences (the representation in terms of the Bessel functions already has 
Minkowski space contributions subtracted). The 5-function shell generates a jump in the 
normal derivative of the eigenfunctions; call the jump U (here, U = A/a). The leading 
heat-kernel coefficients for this situation, namely for functions which are continuous across 
the boundary but which have a jump of the first normal derivative at the boundary, have 
been determined in Ref. |36|; the relevant coefficient is given in Theorem 7.1, p. 139 of 
that reference. The results there are very general; for our purpose there is exactly one term 
that survives, namely 

a 4 = — [ dlU 3 (7.6) 
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which shows that 

So no matter along which line the 5-function shell is concentrated, the first two orders in 
a weak-coupling expansion do not contribute any divergences in the total energy. But the 
third order does, and the divergence is given above. 

For the example considered, as mentioned, U = A/a is constant, and the integration leads 
to the length of the line which is 2ira. Thus we get for this particular example 

«=9^ < 7 ' 8) 

This can be easily checked from the explicit representation we have for the energy. We have 
already seen that the first two orders in A identically vanish, while the part of the third 
order that potentially contributes a divergent piece is 

1 00 poo J 1 

— J2 / dxx 2 - 2s --X 3 KUx)I 3 M. (7.9) 
rcr ^— ' J ax 3 



m=—oo 



The m = contribution is well behaved about s = 0; while for the remaining sum using 

Kl{mz)ll{mz) ~ 1 * /2 , (7.10) 

8m- 3 (1 + z z J ,i ' z 

we see that the leading contribution is from f)6.7j) 

^LVm 2 ~ 2s I" dzz 2 - 2s — — - 

Una 2 ^ J dz8m 3 (l + z 2 ) 5 / 2 

m=l u v ' 

A 3 „ ^ _ , f°° , 2 __ 2s d 1 



dzz 



9QTia 2SUK ' J dz(l + z 2 yi 2 

* fa(1 + 2 /P-f + » * +0(A ( , n) 



967ra 2 ^ v 7 T(3/2) 96vra 2 s 
in perfect agreement with the heat-kernel prediction ()7.8|) . 



VIII. STRONG COUPLING 

The strong-coupling limit of the energy (|3.5|) . that is, the Casimir energy of a Dirichlet 
cylinder, 



^ D = -^ £ / dxx 2 —\nl m (x)K n 
Sira z ^ J dx 



(8.1) 



m=— 00 



was worked out to high accuracy by Gosdzinsky and Romeo [31| 



0.000614794033 



It was later redone with less accuracy by Nesterenko and Pirozhenko j3JJ. (Note that the 
preprint version of the latter contains several internal sign errors.) 
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For completeness, let us sketch the evaluation here. We carry out a numerical calculation 
(very similar to that of in the spirit of Sec. IV Bl We add and subtract the leading 
uniform asymptotic expansion (for m = the asymptotic behavior) as follows: 



£ 



D 



di^jc i oc 



ln(2xI (x)K (x)) - - ' 



1 + x 2 



m=l 



dxx 2 — 
dx 



f xt\ 1 ri(t) 



\n(2xI m (x)K m (x))-\n[ — )- - 



\ m 



2 m 2 



- 2 



I OO \ poo I 00 PC 

m=lj J{> m=l JU 



d 



dxx — In xt 

dx 



°° 2 d 
doc oc ■ 

dx 



m=l 

1 (\ ^ 



n(t) 



m=l/ 



dxx z 



Al + x 2 

d 1 

cfe 1 + X 2 



(8.3) 



In the first two terms we have subtracted the leading asymptotic behavior so the resulting 
integrals are convergent. Those terms are restored in the fourth, fifth, and sixth terms. 
The most divergent part of the Bessel functions are removed by the insertion of 2x in the 
corresponding integral, and its removal in the third term. (Elsewhere, such terms have 
been referred to as "contact terms.") The terms involving Bessel functions are evaluated 
numerically, where it is observed that the asymptotic value of the summand (for large m) 
in the second term is l/32m 2 . The fourth term is evaluated by writing it as 



21im V s 



m 



,2-s 



m=l 



As 



dz- 



1 + Z 2 



2C'(-2) 



03) 
2tt 2 



(8.4) 



while the same argument, as anticipated, shows that the third "contact" term is zero, 1 while 
the sixth term is 



— lim 

2 s^o 



CM + 5 



1 1 , 

- = -ln27T. 

s 4 



■5) 



The fifth term is elementary. The result then is 
1 



£ 



D 



Ana 2 
1 



(0.010963 - 0.0227032 + + 0.0304485 + 0.21875 - 0.229735) 
(0.007724) = 0006146 



Ana 2 a A 
which agrees with ()8.2|) to the fourth significant figure. 



(8.6) 



This argument is a bit suspect, since the analytic continuation that defines the integrals has no common 
region of existence. Thus the argument in the following subsection may be preferable. 
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A. Exponential Regulator 



As in Sec. El it may seem more satisfactory to insert an exponential regulator rather 
than use analytic regularization. Now it is the third, fourth, and sixth terms in ()8.3j) that 
must be treated. The latter is just the negative of (|5.14j) . We can combine the third and 
fourth terms to give using (|5.13|) 



1 



2 



dz z 3 d 2 



1 



z 2 + 5 2 dz 2 e z — 1 ^ 

The latter integral may be evaluated by writing it as an integral along the entire z axis, and 
closing the contour in the upper half plane, thereby encircling the poles at i5 and at 2inn, 
where n is a positive integer. The residue theorem then gives for that integral 

so once again we obtain the same finite part as in ()8.4jl . In this way of proceeding, then, in 
addition to the finite part in (J8.6)) . we obtain divergent terms 

111 . , 

(8.9) 



c div 



+ 



64a 2 <5 8na 2 5 2 4a 2 <5 3 ' 
which, with the previous interpretation for 8, implies divergent terms in the energy pro- 
portional to L/a (shape), L (length), and ah (area), respectively. Such terms presumably 
are to be subsumed in a renormalization of parameters in the model. Had a logarithmic 
divergence occurred [as does occur in weak coupling in 0(A 3 )] such a renormalization would 
apparently be impossible — however, see Sec. |Xj 



IX. LOCAL ENERGY DENSITY 



We compute the energy density from the stress tensor (|3.1j) . or 



(T 00 ) = ^(9°9 0/ + V- V')G(x,x') 

— t 



-V 2 G(x,x) 



i rev r 00 °° [ / rn2 \ 



m=— oo 



2£-d r rd r g(r, r) 



(9.1) 



We omit the free part of the Green's function, since that corresponds to the energy that 
would be present in the vacuum in the absence of the cylinder. When we insert the remainder 
of the Green's function (|2.10bj) . we obtain the following expression for the energy density 
outside the cylindrical shell: 



\ POO POO °° 

u(r) = (T™ -T$) = - j^ J d( J dk j 



IKko) 



m=—oo 



X 



m 



+ XI m (Ka)K m (Ka) 
1 d d 



2oo 2 + k 2 + —\ K 2 m {Kr) + K 2 K£(Kr) - 2^-—r—K^r 



r dr dr 



r > a. 



(9.2) 
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The last factor in square brackets can be easily seen to be, from the modified Bessel equation, 

2^( K r) + l^ill.r§-Kl( K r). (9.3) 

For the interior region, r < a, we have the corresponding expression for the energy density 
with I m <-> K m . 



A. Total and Surface Energy 



We first need to verify that we recover the expression for the energy found in Sec. IIVI So 
let us integrate expression ()9.2j) over the region exterior of the cylinder, and the correspond- 
ing interior expression over the inside region. The second term in (jfj.Hjl is a total derivative, 
while the first may be integrated according to (j4.2b|) . In fact that term is exactly the one 
evaluated in Sec. IIVI The result is 



/■OO POO J 

2tt / drru(r) = V / dxx 2 — In [1 + \I m {x)K m {x)\ 

JIJ m=-oo J ° 

poo 00 

- (1 - 4 «»wi dxx E I 



I m (x)K m (x) 



„ ; __ x +XI m (x)K m (x) 

(9.4) 

The first term is the total energy (|3.5|) . but what do we make of the second term? In 
strong coupling, it would represent a constant that should have no physical significance (a 
contact term — it is independent of a if we revert to the physical variable k as the integration 
variable). In general, however, there is another contribution to the total e nergy resi ding 
precisely on the singular surface. This surface energy is given in general by E3, 13, H H 



£ = -^r-^ / dS ■ VG{x, x') 
2« Js 



(9.5) 



where the normal to the surface is out of the region in question. In this case it is easy to see 
that (£ exactly equals the negative of the second term in ()9.4|) . This is an example of the 
general theorem 

(dr)w(r) + £ = E, (9.6) 



that is, the total energy E is the sum of the integrated local energy density and the surface 
energy. A consequence of this theorem is that the total energy, unlike the local energy 
density, is independent of the conformal parameter £. For more on the surface term, see 
Sec. IX Bl and Appendix B. 



B. Surface Divergences 

We now turn to an examination of the behavior of the local energy density (J9.2)) as r 
approaches a from outside the cylinder. To do this we use the uniform asymptotic expansion 
()5.7a|) . ()5.7b|) . where now we need to know that drj/dz = 1/zt. Let us begin by considering 
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the strong-coupling limit, a Dirichlet cylinder. If we stop with only the leading asymptotic 
behavior, we obtain the expression 



u(r) 



l ,00 - 

o / dn k 2 > 1 

8tt 3 J ^ 

u m=l 



where 



+ (l-40« 2 



A 



7T 1 

2mt z 2 



-k 2 + (1 - 40 k 2 + 



(A -> 00), 



7Tt 

2m 



77(2) — i] (z— 



r / J 



(9.7) 



(9.8) 



and we have carried out the angular integral as in (|4.4j) . Here we ignore the difference 
between r and a except in the exponent, and we now replace k by mz/a. Close to the 
surface, 

2 r — a 



A ~ 7 



i r 

and we carry out the sum over m according to 

00 ^ 1 

2Vm 3 e- m *~-2- 
z / 

m=l 

Then the energy density behaves, as r — > a+, 

3 



r - a < r, 



12 

dx 3 X A 4 



3 Ft 



4^,4 



4 (r - a) 4 ' 



(9.9) 



(9.10) 



u(r) ~ 



647T 2 (r 
1 1 



16tt 2 (r - a) 4 



-— / rfz^ 5 + t 3 (l-80] 
1-60- 



(9.11) 



This is the universal surface divergence first discovered by Deutsch and Candelas \2L It 
therefore occurs, with precisely the same numerical coefficient, near a Dirichlet plate |32j or 
a Dirichlet sphere 11]. It is utterly without physical significance, and may be eliminated 
with the conformal choice for the parameter £, £ = 1/6. 

We will henceforth make this conformal choice. Then the leading divergence depends 
upon the curvature. This was also worked out by Deutsch and Candelas [2]; for the case of 
a cylinder, that result is 



u(r) 



7207T 2 r(r — a) J 



a+, 



(9.12) 



exactly 1/2 that for a Dirichlet sphere of radius a [111 ]. Here, this result may be straightfor- 
wardly derived by keeping the 1/m corrections in the uniform asymptotic expansion (|5.7ajh 
()5.7b[) . as well as the next term in the expansion of Xi 



2r — a I r — a 
+ t 



X 

t r \ r 

(Note that there is a sign error in (4.8) in Ref. |llj.) 



(9.13) 



19 



C. Weak Coupling 

Let us now expand the energy density ()9.2j) for small coupling, 

1D7T J — J -co „_ n 



X 



m=— 00 ?i=0 
,2 



-« 2 + (1 - 40 /t 2 + 



^f) Kl( K r) + (l-^ 2 K'^r)}. (9.11) 



If we again use the leading uniform asymptotic expansions for the Bessel functions, we obtain 
the expression for the leading behavior of the term of order A n , 



u {n) (r) ~ 



4-j ~ / dzzY^m^e-^r- 1 ^ + 1-80- (9.15) 
Tr V l J Jo TO=1 

The sum on m is asymptotic to 

00 / + \ 4— n 

^ m 3-n e - mX ^ (3 _ n)! (_Z_) ; r ^ fl+) (Q16) 



m=l 



so the most singular behavior of the order X n term is, as r — > a+, 

,<»)^ _ , ^»(3-n)!(l-6e) 



uW(r) ~ (-A) n v - V v (9.17) 



This is exactly the result found for the weak-coupling limit for a 5-sphere jUj and for a 
5-plane j34T ]. so this is also a universal result, without physical significance. It may be made 
to vanish by choosing the conformal value £ = 1/6. 

With this conformal choice, once again we must expand to higher order. Besides the 
corrections noted above, in (|5.7a|) . (|5.7b|) . and (|9.13jl . we also need 

t = tiza/r) ~ t + (t-t 3 ) 1 —^, r^a, (9.18) 

r 

Then a quite simple calculation gives 

» W ~ r^a + , (9.19) 



28807r 2 r n+1 ('r - a 



which is analytically continued from the region 1 < Ken < 3. Remarkably, this is exactly 
one-half the result found in the same weak-coupling expansion for the leading conformal 



divergence outside a sphere |ll(. Therefore, like the strong-coupling result (|9.12|) . this limit 



is universal, depending on the sum of the principal curvatures of the interface. 



X. CYLINDRICAL SHELL OF FINITE THICKNESS 

In this section we regard the shell (annulus) to have a finite thickness 5. We consider the 
potential 

Ant = ~^Mr), (10.1) 
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where 



0, r < a_, 
rr[ r) . //. a < i- < (i . . 
0, a i < r. 



(10.2) 



Here a± = a ± 8/2, and we set M = 1. In the limit as 5 — > we recover the 5-function 
potential considered in the rest of this paper. 

As in Ref. jllj it is straightforward to find the Green's function for this potential. In 
fact, the result may be obtained from the reduced Green's function given there by an evident 
substitution. Here, we content ourselves by stating the result for the Green's function in the 
region of the annulus, a_ < r, r' < a + : 

gm(r,r') = J m (fi;V < )ir m (fi;V > ) + AJ m (K / r)/ m (/tV / ) 

+ B[I m {K'r)K m {K'r') + A m («V)/ m («V)] + CK m ^'r)K m (K'r'), (10.3) 

where k 1 = \J 'k 2 + Xh/a. The coefficients appearing here are 

1 



A = ——[Kl' m (Ka-)K m (K'a-) — n'I m (K,a^)K' m (K,'a)] 

y.[KK' m (Ka + )K m (K a + ) - K,'K m (na + )K' m (Ka + )], 

B = —[Kl' m (Ka )I m (KaJ) — k I m (Ka_)I' m (K a )} 

x[KK' m (Ka + )K m (K'a + ) - K'K m (Ka + )K' m (K'a + )), 

C = — — [Kl' m (KaJ)I m (K'a-) — K'I m (Ka_)I' m (n'a_)} 

x[KK' m (Ka + )I m (K'a + ) - n'K m (Ka + )I' m (K,'a + )], 



(10.4a) 
(10.4b) 
(10.4c) 



where the denominator is 



[Kl' m (Ka_)K m (n' a_) - K'I m (KaJ)K' m (K'aJ)] 
x[KK' m (Ka + )I m (K'a + ) - K'K m (Ka + )I' m (K f a + )] 
- [Kl' m (Ka-)I m (ri a_) - K'I m {Ka-)l' m {K'aJj\ 
x[KK' m (Ka + )K m (K'a + ) - «' 'K m (Ka+)K' m (rfa+)\. 



(10.5) 



The general expression for the energy density within the shell is given in terms of these 
coefficients by 



u 



dKK ? 



k m , 9 Xh , . , 9 

+ — + k 2 + 4£k' 2 

2 r z a 



m 



n' 2 r 2 



x[AI 2 m (n'r) + CKlXn'r) + 2BK m (K'r)I m (^r)} 

+ - 40K(«/r) + CK*{k't) + 2BK> m {K'r)I' m {K'r)\ 



d,K K 



j | 1-4£ 1 d^d 
2 r dr dr 



x 



Wl^'r) + CKl(K'r) + 2BK m (K l r)I m (K l r)}. 



(10.6) 
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A. Leading Surface Divergence 



The above expressions are somewhat formidable. Therefore, to isolate the most divergent 
structure, we replace the Bessel functions by the leading uniform asymptotic behavior (|5.7j) . 
A simple calculation implies 



A 
B 

C 



U - t' 



-2mr)', 



t + +t' + 

t + -t' + t_- t'_ 
t+ + f + tl + t'_ 
t- ~ t'_ 



2m(r)'_-v'+) 



t_ + t'_ 



2mrf_ 



(10.7a) 
(10.7b) 
(10.7c) 



using the notation t + = t(z + ), rj'_ = T)(z'_), etc., where for example z'_ = n'd-fm. If we now 
insert this approximation into the form po.fi j) for the energy density, we find 



u 



(T uu ) 



i ^ r 

+ m=l J{> 



X 



t + + t' + 

~m 2 z 2 



1 - 80 + 



t- + 1'_ 

Xha 2 , 
± + 



m 2 a 2 l _ 



;i-4o 



m 2 zl U $- % e 2m «- "<) 



't+ + t', t_ + f_ 



(10.8) 



If we are interested in the surface divergence as r approaches the outer radius a + , the 
dominant term comes from the first exponential factor only. Because we are considering the 
limit Xha <C m 2 , we have 



t' 



f | 1 Xk a l f 2 

1+11 2^T t+ '' 



and we have 



u 



Xh/a 
327r 2 a 2 



oo /*oo 

J2 m dzzt(l -8^ + t 2 )e 2mW -"V). 

m=l ^° 



(10.9) 

(10.10) 



The sum over m is carried out according to (|9.1fij) . or 



OO / 

J2 me 2m ^-^ ~ ( 

7Tl=l ^ 



rtL 



2(r - a. 



(10.11) 



and the remaining integrals over z are elementary. The result is 

Xh 1 - 6f 



967r 2 a (r — a + ) 2 



r->ai 



10.12) 



This is the expected universal divergence of a scalar field near a surface of discontinuity 
which is without significance, and which may once again be eliminated by setting £ = 1/6. 
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B. Surface Energy 



In this subsection we want to establish that the surface energy (E ()9.5|) is the same as the 
integrated local energy density in the shell when the limit 5 — > is taken. To examine this 
limit, we consider Xh/a 3> k 2 . So we apply the uniform asymptotic expansion for the Bessel 
functions of k' only. We must keep the first two terms in powers of k <C k': 



, 2 I m (Ka-)K m (Ka + ) . , , 

-k : — : — , — sinn m{T]_ — rj + 



mz 



K K 

m 



1 t 1 



1 t' 



cosh m(i]'_ — 



Because we are now regarding the shell as very thin, 

5 1 



where 



using the Wronskian ()2.7|) we get 



a t 



rn 



/ • 



VXhc 



- — [1 + XI m (Ka)K rn (Ka)}. 
a 



Then we immediately find the interior coefficients: 

A 

1 



I m {Ka)K m (Ka) 



~ —y/Xha ; : 

2 1 + XI m (Ka)K m (Ka) 



-2mr)' 



I m (Ka)K m (K,a) 



B ~ -VXha — — j — — — j — r, 

2 1 + XI m {Ka)K m {Ka) 

1 r— — I m (Ka)K m (Ka) 

C ~ —VXna- — — — — - 

Zir 1 + XlmiKajKjniKa) 



(10.13) 

(10.14) 
(10.15) 

(10.16) 

(10.17a) 
(10.17b) 
(10.17c) 



We now insert this in (jK).fij) and keep only the largest terms, thereby neglecting k 2 relative 
to Xh/a. This gives a leading term proportional to h, which when multiplied by the area of 
the annulus lisab gives for the energy in the shell 



£ ann = 2-nbau ~ (1 - 4£) 



A 



4ira 2 



E 



m=— oo 



dna Ka- 



I m (Ka)K m (Ka) 
1 + XI m {Ka)K m (na) 



10.18) 



which is exactly the form of the surface energy given by the negative of the second term in 

(HE1. 

In particular, note that the term in € of order A 3 is, for the conformal value £ = 1/6, 
exactly equal to that term in the total energy £ (|3.5|) : [see 1)7.9)1 ] 

€ (3) = £ (3) . (10.19) 

This means that the divergence encountered in the global energy (J7.8j) is exactly accounted 
for by the divergence in the surface energy, which would seem to provide strong evidence in 
favor of the renormalizablity of that divergence. 
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XI. CONCLUSIONS 



In this paper we have extended the considerations applied to a spherical geometry in 



Ref. 1 1] to a cylindrical one. Results for the general structure of the local and global energies 
are rather as expected, with the leading conformal divergences in the local energy density as 
the surface is approached being reduced from their spherical values by a factor of 1/2. The 
only possibly surprising result is that the weak-coupling limit of the total Casimir energy for a 
5-function cylindrical shell vanishes through order A 2 , in agreement with other perturbative 
results for a c ylin der, as seen in Table |U We do not yet have a complete explanation of 
this. Schaden |39| has presented a semiclassical calculation in terms of periodic rays which 
sheds light of this general phenomenon. His technique closely reproduces Boyer's result 
for the spherical shell, while giving a vanishing result to all orders for the cylinder. This 
cannot be the whole story, however, since Casimir energies for cylinders only vanish through 
second order, and do not vanish for a perfectly conducting cylinder. Here, the weak-coupling 
Casimir energy diverges in 0(X 3 ), which has its origin in the singular nature of the surface 
energy. 

APPENDIX A: PERTURBATIVE DIVERGENCES 

A conventional field theorist might be surprised that we are able to extract finite (actually 
zero) values for the global energy in order A and A 2 , when the latter contributions would 
seem to correspond to divergent Feynman graphs. Thus the first-order energy is given by 
(T is the very large time interval characterizing the space-time volume under consideration) 

= ~D + (Q) / (dx)a(x), (Al) 

which involves the massless scalar propagator D + (x) at zero coordinate, corresponding to 
equal field points. The latter is ill-defined. For example, for our cylindrical 5-shell geometry, 
we would compute from this the energy per unit length 

£« = £i e ->0+, (A2) 

which, although divergent, is imaginary and independent of a. Such a contribution should 
be irrelevant. 

Less divergent, but still ambiguous, is the second-order expression 

E {2) = ^ J(dx)(dy)a(x)D + (x-y)a(y)D + (y-x). (A3) 



This may be evaluated by a method very similar to that employed in Ref. [32]. If we work 
in D + 1 dimensions, we obtain 

\2 P ( 3—1) ^ -r poo poo 

^—T ^o-LJ^ l «** D - 2 Jt(*)[ *uu*-w (l - ur -w. (A4) 

For D > 1 the second integral may be evaluated as a beta function. The integral over the 
squared Bessel function is, for 1 < Re.D < 2, 

2 r i 1 ~ ~ ) r i^r) 

dxx J (x) = - /3-D\ • 

^V 71 " 1 \—2~) 
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When this is inserted into (|A4|) a factor of T ((3 — D)/2) remains in the denominator, so the 
second-order energy evidently vanishes in three dimensions. 

We can also proceed directly in three dimensions, if we regulate the integral using a 
proper time method as described in Ref. 32||. If, as there, we exclude a coincident field term 
(contact term), we are led to 



£(2) 



A 2 7T 2 d 

4 da 



POO 

/ dqq 2a+1 J 2 (qa) 
Jo 



(A6) 



a=0 



We encounter the same Bessel-function integral as in (jA5|) . so continuing the integral from 
— 1 < Re a < —1/2 we again see that £^ = 0. (This same continuation process gives 
the correct second-order energy for a sphere: = \ 2 /32ira for a dimensionless coupling 
constant A.) 

In the balance of this appendix we sketch an evaluation of along the lines of that 
given for 8^ in Sec. IV Bl We regulate the sum and integral by inserting an exponential 
cutoff, 5^0+: 



A 



Ana? 

which we break up into five parts 



^ + [ dxx 2 -^I m (x)K m (x)e 

771=1/ 



T + II + III + IV + V). 



(A7) 



8vra 2 



(A8) 



The first term is the m = contribution, suitably subtracted to make it convergent (so the 
convergence factor may be omitted), 



I 



d 

dx 



I (x)K()(x 



2VT+ 



-1. 



The second term is the above subtraction, 



n = i 

2 



d 



dx VI + x 2 



e ~ 



-26 + 1 ' 



(A9) 



(A10) 



as may be verified by breaking the integral in two parts at A, 1 A <ti 1/5. The third term 
is the sum over the mth Bessel function with the two leading asymptotic approximants ()6.3|) 
subtracted: 



rn=l 



x 2 d 



dxx' 



dx 



I m (x)K m (x) - 



t 

2m 



1 



5m z 



'1 - 6t 2 + 5t 4 ' 



0. 



(All) 



Numerically, each term in the sum seems to be zero to machine accuracy. This is verified by 
computing the higher-order terms in that expansion, in terms of the polynomials in (|5.10J) . 



I m (x)K m (x) - — 
2m 



1 



5m z 



4m 5 



r 2 (*) - ^r 2 (t) 



4m 7 



'l-6t 2 + 5r 



+ 



(A12) 
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which terms are easily seen to integrate to zero. The fourth term is the leading subtraction 
which appeared in the third term: 

iv = J2 m I dzz2 ij-A e ~ mzS - ( A13 ) 

m=l V / 

If we first carry out the sum on m we obtain 

1 f°° 1 1 

W = ~ dzz ? 



4 J (1 + z 2 fl 2 sinh 2 z5/2 

~-^ + T6-l (A14 » 

as again may be easily verified by breaking up the integral. The final term, if unregulated, 
is the form of infinity times zero: 



8 ^ m In 

m=l J{i 



oo j 

d z z 2 jAf - 6t 5 + 5t 7 )e- mzS . (A15) 



Here the sum on m gives 



oo 

J2~ e ~ mZ 5 = -ln(l-e^), (A16) 



m 

m=l 



and so we can write 



V=-| f&fl- uru' 2 -*^ 2 - Qu b ' 2 + 5u 7 ' 2 ) 
16 da J 

Adding together (JMJ), (fATojl . f[XTT)l . (fA"T4jl . and (fATTjl . we obtain 



a=0 



1 

6' 



(A17) 



^ (1) = ^A^ + ' ( A18 ) 

that is, the 1/5 and constant terms cancel. The remaining divergence may be interpreted as 
an irrelevant constant, since 5 = r/a, r being regarded as a point-splitting parameter. The 
correspondence of the terms in this evaluation with that Sec. IVII is rather immediate. 



APPENDIX B: SURFACE ENERGY 



Here, we suggest an alternative derivation of the surface energy term for £ = 0. For the 
scalar Lagrangian 

C = <j>dy - ±5(r - a)0 2 , (Bl) 

the response to a general coordinate transformation yields the expected stress tensor, in- 
cluding the interaction term with the cylindrical surface: 

T"" = d»(f>d v (f> + g^C. (B2) 
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The integral of the ^-function term in the energy density should be the surface energy 

f°° A 

£ = 2n drr— <5(r - a)(0 2 (r)> = 7rA(0(a) 2 ), (B3) 
Jo 2a 

or, in terms of the Green's function 

VTA f°° du f°° dk 1 ^ . . 



X [""dxx T X JSMm (B4) 
Jo S-l + \I m {x)K m {xY 1 ; 



Ana 2 



m=—oo 



This is similar to the surface energy given by (jIU.ISp for £ = 0; the difference between the 
two expressions is 



A 



Aira 2 

m=—oo 



/ dxx I m (x)K m (x) = -— dKKK (0) 



(B5) 



if we use the addition theorem (|5.2jl in the singular limit p' — > p. Such a term, a singular 
constant, would seem to be a physically irrelevant contact term, so the two versions of the 
surface energy appear to be equivalent. 
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